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ABSTRACT 

Aims. We derive the dense core structure and the water abundance in four massive star-forming regions which may help understand the 

earliest stages of massive star formation. 

Methods. We present Herschel-HIFI observations of the para-H20 in - Ooo and 2o2 - 1 1 1 and the para-H2**0 in - Ooo transitions. The envelope 

contribution to the line profiles is separated from contributions by outflows and foreground clouds. The envelope contribution is modeled using 

Monte-Carlo radiative transfer codes for dust and molecular lines (MC3D and RATRAN), with the water abundance and the turbulent velocity 

width as free parameters. 

Results. While the outflows are mostly seen in emission in high-J lines, envelopes are seen in absorption in ground-state lines, which are 

almost saturated. The derived water abundances range from 5x10'" to 4x10 ** in the outer envelopes. We detect cold clouds surrounding 

the protostar envelope, thanks to the very high quality of the Herschel-HIFI data and the unique ability of water to probe them. Several 

foreground clouds are also detected along the line of sight. 

Conclusions. The low H2O abundances in massive dense cores are in accordance with the expectation that high densities and low temperatures 

lead to freeze-out of water on dust grains. The spread in abundance values is not clearly linked to physical properties of the sources. 

Key words. ISM: dust, extinction - ISM: molecules - ISM: abundances 
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1. Introduction 

Massive stars (>1OM0) play a major role in the interstel- 
lar energy budget and the shaping of the Galactic environment 
(IZinnecker & Yorka l2007h . However, the formation of such 
high-mass stars is not well understood due to several reasons: 
they are rare, have a short evolution time scale, they are born 
deeply embedded, and are far from the solar system. 

The main sequence lifetime of massive stars is pre- 
ceded by an embedded phase which subdivides into several 
classes of objects: massive pre-stellar cores (mPSC), which 
are loc al temperature minima a nd density maxima within dark 
clouds dSridharan et alJllOOSi) : high-mass protostellar objects 
(HMPO), where a central protostar is surrounded by a massive 
envelo pe with a centrally p eaked temperature and density distri- 
bution dvan der Tak et alJllOOOi) : hot molecular cores (HMC), 
which have larger masses of warm gas and dust, and high abun- 



* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with parti- 
cipation of NASA. 



dances of complex organic molecules which have evaporated 
off dust grains and /or formed by warm gas-phase chemistry 
JMotte et all 120031) : and ultracompact H II regions (UCHII), 
which show large pockets of ionized gas confined to the star 
(Churchwell et al. 1990). A key question is to what extent these 
phases represent differences in luminosity and/or age, and if all 
high-mass stars pass through all these phases. 



The water molecule is thought to be a sensitive tracer of 
physical conditions in star-forming regions, which acts as a nat- 
ural filter for warm gas because of it s large abundance varia - 
tions between hot and cold regions ('van der Tak et al."200^. 
Moreover, because the dust continuum is strong at the higher 
frequencies, water lines connecting with the lowest energy 
levels can be seen in absorption, thus providing an alterna- 
tive metho d to probe different depths in t he protostellar en- 
vironment iPoelman & van de r Taldl2007h . Measurements of 
the abundance of water are therefore a step toward understand- 
ing the energy budget of star-forming regions, and thus of the 
star formation process itself. 
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This paper presents water observations performed with 
the Heterodyne I nstrum ent for the Far-Infrared (HIFI; 
[de Graauw & et a]] 20101) on-boa rd ES A's Herschel Space 
Observatory dPilbratt & et al.ll2010l) . We use the p-HiO ground- 
state line and two lines which constrain the excitation and 
optical depth (Table [1]), all three lying at similar frequencies 
and observed at similar resolution. The sources are four 
massive star-forming regions (Table|2li: the HMCs G3 1.41 +0.31 
and G29.96-0.02 and the HMPOs W33A and W43-MM1. We 
compare our r esults with those for two other regions: the UCHII 
region DR21 jvan der Tak et alj |20Tol) and the HMPO W3 IRS5 
dChavarria et al.ll20ig, this volume). 

However, the aim is to discover trends in the water line 
emission for future extended studies, identifying links in the 
water abundance between the various evolutionary stages of 
high-mass star formation and to use water as probe of the 
gas dynamics around protostars. Given the small number of 
sources and Unes observed, it is premature to look for general 
trends. The large amount of upcoming Herschel-HIFI data 
will help on this issue. 



Table 1. List of lines. 



Molecule 


Transition 


V (GHz) 


£up (K) 


ttcril'^ (cm -^) 


o-,-„„ (mK) 


H2O 
H2O 


111 - Ooo 
2o2-ln 
111 - Ooo 


1113.343 
987.927 
1101.698 


53.4 
100.8 
53.4 


1.7x10" 
2.1x10" 
1.7x10" 


40 
50 
40 



Notes. '"' Values at 20 K from collision rates of lGrosiean et al.l i2003h 



Table 2. List of sources. 



R.A. 
J2000 



Dec. 
J2000 



L 

(lO-'Lo) 



d" Vlsr 

(kpc) (kms-') 



G3 1. 4 1+0.3 1 
G29.96-0.02 
W33A 
W43-MM1 



I8M7"'34.3" 
18''46"'03.8" 
18''14"'39.P 
18''47'"47.0" 



-01°12'46.0" 
-02° 39' 22.0" 
-17°52'07.0" 
-01°54'28.0" 



15 
20 
8.5 





7.9 +98.8 

7.4 +98.7 
4.0 +37.5 

5.5 +98.8 



Notes. <°' Values from fHatchell & van der Tald (|2003|[), exc ept W43- 
MMl JMotte et al.ll2003h and W33A jvan der Tak et ali200(3h . 



2. Observations 

The four regions have been observed with HIFI on the 3"*, 4* 
and 6* of March 2010 (see Table|2|. Spectra were taken in dou- 
ble sideband mode using receivers 4a (P-H2O at 988 GHz) and 
4b (P-H2O and p-Hi^^O at 1113 GHz and 1102 GHz) with vlo 
= 980 GHz and 1108 GHz respectively. The observations are 
part of the Priority Science Program (PSP) of the Guaranteed- 
Time Key Program Water In Star-forming regions with Herschel 
(WISH; Van Dishoeck et al., in prep.). 

Data were simultaneously taken with the acousto-optical 
Wide-Band Spectrometer (WBS) and the correlator-based High- 
Resolution Spectrometer (HRS), in both horizontal and verti- 
cal polarizations. This paper focuses on data from the WBS, 
which covers 1140 MHz band width at 1.1 M Hz spectral reso- 
lution (~0.3 kms ') faoelfse ma et alJllOlOl) . System temper- 
atures range between 350 K around 1113 GHz and 450 K 



around 988 GHz; receiver 4a in V polarization shows partic- 
ularly high values. Integration times (ON+OFF) were 193 s 
for the 1113 GHz and the 1102 GHz lines, and 206 s for the 
988 GHz line, for each source, and the rms noise levels reached 
are 40-50 mK (Table [TJ. Observations have been reduced with 
the Herschel Interactive Processing Environment (HIPE) ver- 
sion 2.8. The intensity scale is converted to T^^b using main beam 
efficiencies of 0.74. The double-side band continuum level has 
been divided by 2 to make its brightness directly comparable to 
that of the lines, which are measured in single sideband. 

3. Results 

Observed water lines for the four studied regions are shown 
in Fig. [1] The H2O In - Ooo line shows an absorption at 
the systemic velocity (Vlsr) in all sources. In all cases except 
G3 1.4 1+0.31, outflow wings are detected close to the main ab- 
sorption, with a maximal shift of 3 kms ^ These wings are 
seen in emissio n, which indicates an origin i n hot, low-density 
(10^ cm"^) gas dPoelman & van der Takl2007l) . Absorption fea- 
tures are seen over a large velocity range in G29.96-0.02, 
W33A, W43-MM1 and more weakly in G31.41+0.31. The ab- 
sorptions at velocity offsets >4 kms"' likely originate in cold 
foreground clouds on the line of sight to the source. In contrast, 
the absorption features at smaller velocity offsets are plausi- 
bly related to cold clouds surrounding the dense cores (which 
other studies call the protostellar envelopes), which are all part 
of large-scale molecular clouds (see Fig. |2]i. Table |3] presents 
a Gaussian decomposition of the Hne profiles around Vlsr; 
Appendix A shows Gaussian iterative decompositions of 
the absorption profiles of the ground-state transition over 
the full velocity range, showing a large number of veloc- 
ity components thanks to the high velocity resolution of 
the HIFI instrument. The absorptions at Vlsr are saturated 
for G3 1.41 +0.31 and W33A and nearly saturated for the other 
sources, which indicates abundances ar ound -lO*^ for the outer 
cold p arts of the massive dense cores dPoelman & van der TakI 
I2007h . 

The H2O 2o2 - 1 1 1 line always appears in emission and shows 
a broad and a narrower velocity component (Fig.[T]). In addition, 
the spectra of G3 1.41 +0.31 and W43-MM1 show two wefl de- 
fined self-absorption features which appear at the source veloc- 
ity. With its high /iup, this transition mainly traces warm gas, and 
the presence of these absorption features in G3 1.41 +0.31 and 
W43-MM1 suggests a higher water abundance in these sources 
than in G29.96-0.02 and W33A. The components seen in emis- 
sion have Gaussian shapes, one being wider (FWHM=20- 
40 kms^') than the other (FWHM=6.4-8.0 kms^'). We asso- 
ciate the broad component with high-velocity outflows associ- 
ated with the protostar also seen in 1 1 1 - Ooo line emission. 
This component is symmetric with respect to the source veloc- 
ity in G29.96-0.02 and W43-MM1, blueshifted by 2.8 km s"' in 
W33A, and redshifted by 4.4 kms"' in G3 1.41 +0.31 (Table©. 
The narrower (hereafter ' medium' ) component is potentially as- 
sociated with shocked surrounding material where water is re- 
leased in the gas phase. Indeed, shocks occur at the interface 
between jets and the surrounding dense envelope, with a ve- 
locity close to that of the massive dense core. Similar re sults 
are found in lKristense n et all d2010t) . [Johnstone et al.l d2010h and 
iChayarria et al. (2010). 

The H'^^O 111 - Ooo transition is seen in absorption at the 
source velocity in G31.41 and W43-MM1, which is not sat- 
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Fig.l. Herschel/HIFI spectra of the H2O In - Oqo (top), H2O 2o2 - In (middle) and H^^^O In - Ooo (bottom) lines. Dashed lines 
drawn at Vlsr- 

Table 3. Gaussian decomposition of the line profiles at velocities close to Vlsr- 



Source 


para- 


HjOdi, 


-Ooo) 


para 


H2O(202 


-111) 


para-H^»0(lii 


-Ooo) 




Vlsr 
(kms"') 


Tmb 

(K) 


AV 
(kms"') 


Klsr 
(kms"') 


Tmb 

(K) 


AV 
(kms"') 


Vlsr T^b 
(kms"') (mK) 


AV 
(kms"') 


G31.41+0.31 


95.1 


0.94* 


3.7 


94.6 
99.3 
103.7 


1.37 

0.42* 

0.2 


6.4 
14.0 
-40 


99.5 0.27* 


5.2 


G29.96-0.02 


91.3 
98.5 
99.4 
103.2 


0.26* 
0.90 
0.99* 
0.40* 


3.9 
18.8 
8.4 
2.3 


97.8 
98.2 


1.10 

3.21 


21 
8.0 


98.5 290 


6.0 


W33A 


35.9 
43.0 


0.85* 
0.53 


11.0 
20.0 


35.5 
38.3 


0.34 
1.87 


28 
7.0 


34.2 0.23* 

37.3 280 


11.8 
8.4 


W43-MM1 


98.7 
103.3 


0.87* 
0.43 


13.5 
14.0 


99.6 
99.7 


0.89 
0.31* 


22 
6.8 


99.4 0.19* 


8.7 



Notes. '** Absorption lines are indicated in T.^],JT^ay,tinuum scales. 



urated. This feature originates from the massive envelope. In 
G29.96, the pure emission profile of this transition implies a 
warm diffuse gas origin. Since G29.96 is also the brightest 
source in the narrower component of the 988 GHz line, we sug- 
gest that the Hj^O emission is dominated by shocks at the in- 
terface between jets and the envelope. The P-Cygni-like profile 
for W33A is a mix of both behaviors: the sum of an absorption 
feature due to the massive core and emission from an outflow. 
The fact that the emission is seen only on the red-shifted part 
of the profile is consistent with the outflow components seen in 
1113 GHz and 988 GHz lines, which are also more powerful on 
the red part of the spectrum. 



4. Discussion and conclusions 



To derive the water abundance in the four massive dense cores 
we have removed features related to outflows and foreground 
clouds from the spectrum before any line modeling. The high 
spectral resolution of HlFl is essential in this process, in partic- 
ular for the absorbers with velocities close to that of the central 
source. Studying the Hl,^0 In - Ooo transition prior to the others 
also facilitates to disentangle the envelope contribution, since 
this line is not saturated having a lower optical depth than 
the main H2O isotope. 



M.G. Marseille et al.: Water abundances in high-mass protostellar envelopes: Herschel observations with HIFI 



Once the main contribution is extra cted, we model its pro - 
file according to the method described in lMarseille et alJ (1200 8h : 
first, the dust emission from the massive de nse core is repro - 
duced with the MC3D radiative transfer code (IWolf et al.lll999h . 
including total luminosity and density profile from the literature 
(power-law index p - -1.5); second, the temperature profile 
obtained is used t o model th e line emission with the RATRAN 
code (Hog erheiide & van der Tak 2000) . The free parameters 
are: Xh,o the molecular abundance relative to H2, and ytuib, the 
turbulent velocity width. 

Good fits are obtained for the H^^O 1 1 1 -Ooo transition, which 
is not saturated unlike the H2O lines. The fitting considers both 
the line strength (area and width) and the profile shapes. We have 
computed a grid of Zh,o and Wturb values, adapting step by step 
the grid around the best y^- Using a '^O/'^O ratio of 500, we 
proceed to model the main isotopic water lines. The H2O abun- 
dance is kept constant in our models. We have tried models with 
an abundance increase in the inner region where T > 100 K, but 
the current data do not favor these models above the constant- 
abundance models. 

We estimate the absolute uncertainty in the retrieved H2O 
abundance to be a factor of 10. Since we use the same 
mo deling strategy as the studies by Ivan der Tak et al.l (1201 Ol) 
and IChavarria et al.l ( 1201 Ol) . the abundances obtained should 
be comparable to better than a factor of 3. Our observed 
spread in abundances of a factor of ~100 is much larger than 
this uncertainty. The same range of abundances is found in 
other HIFI-based studies of high-mass star -forming regions 
(Ivan der Tak et al.ll2010l: Chavarria et al.l 20101) . and also in pre- 
vious work w ith ISO ( Boonman et alj2003h and from the ground 
dvan der Tak et al. 200^1 

In conclusion, for the massive star forming regions de- 
scribed in this letter, we clearly detect the contribution of 
the envelope within the dense core. It is limited to a strong 
self-absorbed feature, mainly seen in the ground-state Une. 
In order to evaluate it, we first have to remove emission from 
outflow shocks and absorption by foreground clouds along 
the line of sight. The velocities of the absorbers indicate that 
some are part of the close environment of the source, whUe 
others are physically unrelated. The derived massive dense 
core abundances suggest a strong freeze-out of water on dust 
grains, and imply that water plays only a minor role in the 
thermal balance of the gas. 

The H2O line profiles do not seem to depend on the sup- 
posed evolutionary stage of the source. For example, the two 
'hot molecular cores' G31.41 and G29.96 show very different 
line profiles, and also their H2O abundances differ by a factor 
of ~100. Also, the abundance variations that we have found do 
not seem related to the luminosity of the sources, their temper- 
ature or their turbulent velocity field. However, the number of 
cases treated is not sufficient for a statistical treatment. Future 
studies following the same procedure with a larger number of 
sources should conclude on this issue. Within our sample, the 
highest H2O abundances are derived for G31.41 and W43-MM1 
which show self-absorbed 2o2- 1 11 line profiles (Fig.[T]). As these 
sources are not the most luminous, hot, or active ones in our sam- 
ple, the origin of such a high abundance is unclear 

Firm conclusions about a link between water emission 
behaviour and the evolutionary stage of the source are lim- 
ited by the small number of sources. Our data show that 
water is a useful tool to understand the gas dynamics in 
and around massive stars forming regions. Future multi-fine 
studies of larger samples are highly promising to answer key 
questions about the formation of massive stars. 
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Fig. 2. Extraction of the saturated absorption of para-H20 1 n -Ooo line in W43-MM1. Original profile appears in black bold, residual 
in green bold. 

Table 4. Model parameters and derived water abundances. 



Source 



(Mo) 



Mnin 

(AU) 



Mnax 

(AU) 



"('"mm) "(^nax) 



(cm-3) 



(cm-3) 



(K) 



(K) 



Xu 



"tmb 

(kms-') 



G31.41+0.31 


1500 


200 


22 515 


8.1x10** 


3.1x10" 


406 


43.2 


3.1x10-* 


14 


G29.96-0.02 


700 


200 


20 700 


4.4x10* 


1.9x10* 


489 


50.8 


<5.0xlO-"' 


1.1 


W33A 


4000 


200 


62 000 


3.5x10** 


4.0x10^ 


291 


26.0 


6.0x10-'" 


1.6 


W43-MM1 


2000 


200 


27 500 


5.0x10* 


2.3x10" 


243 


24.7 


4.0x10-* 


3.0 


DR2r 


1650 


2000 


60 520 


1.6x10' 


1.5x10^ 


117 


23.3 


2.0x10-'" 


3.0 


W3-IRS5' 


250 


200 


12 000 


2.9x10* 


2.7x10" 


480 


54.7 


2.0x10-* 


2.0 



Notes. <°^ Values from lvan der Tak et alj r2010h "'' Values from lChavarria et allJIoTol) 



Appendix A: Massive dense core component 
extraction 



The velocity profiles of the H2O 1 1 1 -Oqo line show absorption features at several 
velocities. These absorption features arise in foreground clouds along the line of 
sight or in cold clouds in the neighbourhood of the massive dense core, and are 
not saturated unlike the absorption from the massive envelope. In addition to 
these absorptions, some sources show H2O emission from protostellar outflows. 



This appendix presents our procedure to remove these features in order to 
extract the contribution from the envelope to the line profile. Contrary to others, 
absorption from this part of the object is saturated. We are then able to remove 
other features by iterative Gaussian fits. This process is helped by the high ve- 
locity resolution provided by the Herschel-HIFI instrument, showing accurate 
and "bumpy" profiles in absorptions. Assuming that each bump corresponds to 
a velocity component, they are removed using the Gaussian fitting tool available 
in the HIPE software. Starting from the component with the lowest velocity, they 
are extracted one by one, using the residual of the previous removal to fit the next 
one. This way of fitting insures a very good extraction of velocity component, 
giving a quasi-unique final decomposition of the absoiption features. Results of 
this process are given in Fig.l2l lA~niA.2llA3] and Tables ??, ??, ?? and ??. 
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Fig. A.l. Extraction of the saturated absorption of para-HaO In- Oqo line in W33 A. Original profile appears in black bold, residual 
in green bold. 




Fig. A.2. Extraction of the saturated absorption of para-H20 In- Ooo line in G29.96. Original profile appears in black bold, residual 
in green bold. 
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Fig. A.3. Extraction of the saturated absorption of para-H20 In- Ooo line in G3 1 .4 1 . Original profile appears in black bold, residual 
in green bold. 
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